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Abstract. We study generation of magnetic fields involving large spatial 
scales by time- and space-periodic small-scale parity-invariant flows. The 
anisotropic magnetic eddy diffusivity tensor is calculated by the standard 
procedure involving expansion of magnetic modes and their growth rates in 
power series in the scale ratio. Our simulations, conducted for flows with 
random harmonic composition and exponentially decaying energy spectra, 
demonstrate that enlargement of the spatial scale of magnetic field is benefi- 
cial for generation by time-periodic flows. However, they turn out, in general, 
to be less efficient dynamos, than steady flows. 
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The present work is an extension of the studies carried out by Lanotte et 
al. (2000) and Zheligovsky et al. (2001), who found steady parity-invariant 
flows with a negative magnetic eddy diffusivity to be quite common. A similar 
investigation for turbulent flows is desirable, especially given that flows in 
experimental dynamos are necessarily turbulent (see discussion ibid.). Some 
time-dependent flows are known to be unable to sustain negative magnetic 
eddy diffusivity - namely, the flows, 5-correlated in time at each point in 
space. They advect a mean magnetic field as a passive scalar, and therefore 
eddy diffusivity can only exceed molecular diffusivity (Biferale et al. , 1995). 
Thus the question, whether time-dependent flows can give rise to negative 
magnetic eddy diffusivity, is not trivial. 

However, dynamo simulations for time-dependent turbulent flows are nu- 
merically demanding. For this reason we focus our attention here at an "in- 
termediate" class of flows - those periodic in time, and employ the simplest 
of them: 

v(x,t) = U(x) + v / ^(V c (x)coscjt + V s (x)sincjt). (1) 

Like in the cited papers, a flow is supposed to be 27r-periodic in spatial 
variables x and parity-invariant, i.e. 

v(x,t) = -v(-x,t). (2) 

(Consequently, its space-averaged helicity is zero, and no a-effect is present.) 

Two families of space- and time-periodic flows were closely examined in 
the context of fast kinematic magnetic dynamo theory. Both are general- 
isations of ABC flows, and flows from both families are constructed of a 
small number of spatial Fourier harmonics. A flow of the kind of "modulated 
waves" was considered by Otani (1993) and Childress & Gilbert (1995); it 
belongs to the class (1). "Circularly polarised" flows were employed in simu- 
lations of Galloway & Proctor (1992) and Galloway & O'Brian (1993); near- 
integrable flows of this kind were studied analytically by Ponty et al. (1993, 
1995) and Childress & Gilbert (1995). The assumed time dependence was re- 
sponsible for chaotic behaviour of the flow trajectories (this is necessary for a 
dynamo to be fast), though the flows depended only on two spatial variables 
(for steady flows this rules out chaos). This gave an opportunity to separate 
out dependence on the third spatial coordinate and to make convincing com- 
putations for high magnetic Reynolds numbers R m ~ 10 4 (Galloway & Proc- 
tor, 1992; Galloway & O'Brian, 1993). Circularly polarised non-integrable 
ABC flows depending on three spatial variables were considered by Brummell 
et al. (1999, 2001). 

We consider here a class of flows which seem to be more physically re- 
alistic: like Zheligovsky et al. (2001), we employ flows with a random spa- 
tial harmonic composition and an exponentially decaying energy spectrum. 
(Random-harmonic flows with a slow - hyperbolic - energy spectrum decay 
were also explored ibid., but they were found to be less efficient generators 
of large-scale magnetic field; consequently, such flows are not investigated 
here.) They can be regarded as a "poor man's model" of turbulent flows. 
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Ideally we would like to examine dynamo properties of a large number of 
sample flows from the class under consideration, so that to collect statisti- 
cally sound information. Unfortunately, computations even in the case of 
the simple time dependence (1) are numerically expensive: evaluation of one 
instance of magnetic eddy diffusivity requires more than a day of a CPU of 
a Dec Alpha processor. Hence it appears possible to evaluate only a limited 
number of such instances. However, since spatial Fourier components of the 
flows employed in computations are chosen at random, one can hope that 
results of computations presented here are typical and attributable to a large 
variety of flows. (Note, that flows from the selected class seize to be typical 
if V c = or V s = 0: it is shown in Section 3 that a non-zero contribution 
from the time-periodic part of (1) in the limit of high frequencies requires 
linear independence of the vector fields V c and V s . Thus (1) represents the 
simplest class of time-periodic flows with hopefully a typical behaviour.) 

Also due to numerical complexity of the problem, we have performed 
computations for just one value of molecular magnetic diffusivity, rj — 0.1 . 
The following constraints restrict the choice of 77: On the one hand, if it is 
too large, diffusion dominates, inhibiting generation of magnetic field. On 
the other, if molecular diffusion is smaller than the threshold for the onset of 
generation of a small-scale magnetic field (i.e., magnetic field which has the 
same spatial periodicity, as the flow v), the short-scale instability dominates: 
the large-scale instability becomes "negligible" , since growth rates of large- 
scale magnetic modes are infinitesimally small (large-scale magnetic modes 
being perturbations of small-scale modes with non-vanishing spatial means, 
associated with a zero growth rate). Zheligovsky et al. (2001) explored large- 
scale kinematic dynamos for three values of molecular magnetic diffusivity, 
rj = 0.1, 0.2 and 0.3 (the flow was normalised so that its r.m.s. was 1). For 
all the three values no generation of small-scale magnetic field occurred, and 
for 77 = 0.1 more flows sustained negative magnetic eddy diffusivity, than for 
77 = 0.2 and 0.3 . This has suggested to choose the value rj = 0.1 in our 
simulations; like ibid., we have checked that 77 = 0.1 is above the magnetic 
diffusivity threshold for the onset of generation of a small-scale zero-mean 
magnetic field for every flow, for which magnetic eddy diffusivity has been 
computed. 

In Section 1 we present a mathematical statement of the problem and 
state results of derivation of the magnetic eddy diffusivity tensor (details can 
be found in Appendix). Dependencies of the minimal magnetic diffusivity 
on the ratio of energies of the steady and time-dependent parts of the flow 
and on temporal frequency u are studied numerically. Results of simulations 
are discussed in Section 2. Steady flows are found to be more capable of 
magnetic field generation, than time-periodic ones. Owing to the factor ^/uj 
in (1), in the limit of high temporal frequencies the time-dependent part 
of the flow provides a finite contribution to the magnetic eddy diffusivity 
tensor. In Section 3 the limit magnetic eddy diffusivity tensor is formally 
derived for uj — > 00, and it is shown numerically that flows can have negative 
eddy diffusivity for high frequencies. Our results are briefly summarized in 
the Conclusion. 
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1. Magnetic eddy diffusivity tensor for time-periodic flows 



In this Section we consider a kinematic dynamo problem for a flow v(x, t) 
of time period T, which is 27r-periodic in each Cartesian variable in space, 
solenoidal (V • v = 0) and parity-invariant (2). 

Temporal evolution of a magnetic field h is described by the magnetic 
induction equation 

dh 

— = ??V 2 h + V x (v x h). (3) 

at 

Substituting h = H(x, t)e xt into (3) one finds that a magnetic mode H(x, t) 
satisfies 

AH=-V + W 2 H + Vx(vxH), (4) 
at 

i.e. the kinematic dynamo problem reduces to a Floquet problem for the 
magnetic induction operator. The mode H(x, t) is solenoidal, 

V • H = 0, (5) 

and it is assumed to have the same time period T, as the flow. ReA is then 
the average rate of growth (or decay) of the mode in time. 

We consider magnetic modes, involving large spatial scales. The modes 
are supposed to depend on the fast spatial variable x and on the slow variable 
y = ex. By the chain rule, spatial derivatives in the eigenmode equation (4) 
and the solenoidality condition (5) must be modified: 

V -> V x + eV y (6) 

(the subscripts x and y refer to differentiation in fast and slow variables, 
respectively). The ratio of the two scales, e > 0, is a small parameter of the 
problem. 

Solution to the Floquet problem (4) is sought in the form of power series 

oo 

A = £ (7) 

n=0 

oo 

H = £H„(x,y,t)e" (8) 

n=0 

Substituting (6) and the series (7) and (8) into (4), expanding and equating 
coefficients for each power of e, one obtains a hierarchy of equations. It is 
discussed in Appendix, how all terms of (7) and (8) can be determined by a 
systematic procedure. In particular, it is shown that: 

• The leading term in (7) is A2 and the expansion involves terms with even 
indices only (A = A 2 „+i = for any integer n > 0). 

• All terms of (8) with even indices are parity anti-invariant in fast variables 
(H 2ri (— x, y, t) = H 2n (x, y,t) for any integer n > 0), and all terms with odd 
indices are parity-invariant in fast variables (H 2 n+i(— x, y, t) = — H 2n+ i(x, y, t) 
for any integer n > 0). 



4 



• The leading term in the decomposition of an eigenmode is independent of 
time: H = H (x, y). Its spatial average satisfies the eigenvalue equation for 
the anisotropic magnetic eddy diffusivity operator, 

r?V y (Ho) + V y x £ j2 D ^^^ = ^ (Ho) , (9) 

fc=im=i a y™ 

and the solenoidality condition 

V y • (Ho) = 0. (10) 
Here (•) denotes the mean of a vector field over the cube of periodicity: 

(f) = (2vr)- 3 / f(x,y)dx, 

J[0,2tt]1 

and superscripts enumerate Cartesian components of a vector field. 

Coefficients of the eddy diffusivity tensor can be evaluated using solutions 
of two auxiliary problems: 

" >k + V V 2 S k + V x (v x S fc ) + ^— = (11) 



at oxk 
(the first auxiliary problem), and 



, }j f)]V r m ,i + V x (v x r m , fc ) + 2^ + e m x (v x (S, + e fc )) = (12) 

(the second auxiliary problem); here {e m } is the basis of Cartesian unit vec- 
tors. Then 

1 r T 

D mjfc = — J (v x r m>fc ) dt. 

It can be verified that are parity anti-invariant (Sfc(x, t) = Sk(— x, t)) 
and solenoidal; T m ^ are parity-invariant and satisfy V x • T m ^ + S™ = 0. 

The partial differential operator in the left-hand side of (9) is comprised of 
second-order derivatives with constant coefficients. Consequently, its eigen- 
vectors are Fourier harmonics: (Ho) = he iqy , where q is a (constant) 

wavevector, and h satisfies 

3 3 

qx^E D mife q m h fe + r?|q| 2 h = -A 2 h, (13) 

fc=l m=l 

hq = 0. (14) 

The quantity i] cddy = min| q | =1 (— A 2 ) is regarded as the minimal magnetic 
eddy diffusivity When it is negative, the associated magnetic mode grows 
in time. 
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2. Magnetic eddy diffusivity for flows (1): numerical results 

Numerical simulations have been carried out for flows (1), where U, V c 
and V s are 27r-periodic parity-invariant solenoidal fields. The fields have been 
generated by the procedure which was applied by Zheligovsky et al. (2001): 
(i) a half of Fourier harmonics with random uniformly distributed compo- 
nents are generated (all the rest ones are obtained by complex conjugation, 
so that the vector field is real), (ii) the gradient part of the resultant field is 
projected out, and (Hi) the harmonics are rescaled in each Fourier spherical 
shell to obtain the desirable energy spectrum. The spectrum of flows em- 
ployed in our simulations exponentially decreases by 6 orders of magnitude, 
the Fourier series being truncated at wavenumber 10. The vector fields are 
normalised so that 

£total = 1, (15) 

where 

Eu** = I / V) -3 / 13 M 2 d*ft = / (|U| 2 + i(|V c | 2 + |V S | 2 )) rfx 

1 JO J[0,2n] 3 J[0,2tt] 3 \ 2 J 

is the average total energy of the flow for u = 1 (hence the magnetic Reynolds 
number can be estimated as R m = i]^ 1 ). Equipartition of the average energy 
of the two time-dependent terms in (1) is assumed: 

/ |V c | 2 dx = / |VJ 2 dx. (16) 

J[0,2tt] 3 J[0,2n] 3 

Solutions to the auxiliary problems (11) and (12) are sought in the form 
of Fourier series, with the 64 3 Fourier harmonics resolution in space and 8 
harmonics resolution in time. For this resolution the spatial energy spectra 
of the solutions decay by at least 10 orders of magnitude, and the temporal 
ones - by 4-5 orders of magnitude. 

All computations presented in this paper have been performed for molec- 
ular magnetic diffusivity r] — 0.1 . It has been verified that for each sample 
flow (1), for which a magnetic eddy diffusivity value is reported here, the real 
part of the dominant eigenvalue of the magnetic induction operator acting in 
the space of 27r-periodic small-scale (i.e. independent of the slow variables) 
zero-mean magnetic fields is negative, i.e. r\ = 0.1 is above the magnetic dif- 
fusivity threshold for the onset of generation of a small-scale magnetic field. 
Codes of Zheligovsky (1993a) have been applied for numerical treatment of 
the Floquet problems 3 for small-scale magnetic modes. 

The following questions have been addressed. 

a) How does the minimal magnetic eddy diffusivity r/cddy change when the flow 
(1) is close to a steady one? Magnetic eddy diffusivity has been evaluated for 

3 Note that in the Floquet problem (4) cigenmodes are defined up to a factor e lJuJt , 
where J is an arbitrary integer, the associated eigenvalues differing by iJu>. Hence in 
computations precautions must be taken to converge to the dominant eigenvalue, whose 
imaginary part does not exceed w/2 in absolute value, i.e. whose imaginary part is minimal 
in absolute value in this set of equivalent eigenmodes. 
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temporal frequency u — 1 and molecular viscosity r] — 0.1 for 30 independent 
samples of (1) satisfying (15) and (16), such that the ratio of the average 
energy of the time-dependent part of (1), 



E osc = l -f (|V c | 2 + |V s | 2 )dx, 

2 Jro,27rl3 



'[0: 

to the energy of the steady profile is small: 



'[0,2tt] 3 



|U| 2 dx = 1/400. 



A histogram of the values of 5?7eddy, the amount of change in magnetic eddy 
diffusivity due to introduction of this weak time periodicity into the flow, is 
shown on Fig. 1. Only in 2 cases out of 30 the moderate time dependence 
makes ?7 e ddy to decrease. 



u 



0.00 



0.02 



0.0-1 



0.06 



0.08 



Figure 1. Histogram of values of <5?7eddy, amounts of change in magnetic eddy 
diffusivity due to introduction of weak time periodicity into the flow. 



b) How does minimal magnetic eddy diffusivity depend on how kinetic energy 
is split between the steady and time-dependent parts of the flow? Magnetic 
eddy diffusivity has been evaluated for three sets of sample flows (1) for 
uj — 1. Profiles of the constituent fields U(x), V c (x) and V s (x) are the same 
in each set, and only amplitudes of these fields are varied in such a way that 
(15) and (16) is satisfied. (In fact, in order to compute any of the three 
curves, a sample flow has been chosen out of those used to construct Fig. 1, 
and profiles of its constituent fields have been employed.) 

Graphs of r/eddy versus the ratio E osc /E tota i are plotted on Fig. 2. Though 
graphs representing different flows are quite different in details, Fig. 2 reveals 
a common tendency: a relative increase of energy contained in the time- 
dependent part of the flow is in general accompanied by an overall (although 
not necessarily monotonous) increase of magnetic eddy diffusivity. 
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0.0 0.2 0.4 0.6 0.8 1.0 

Figure 2. Minimal magnetic eddy diffusivity r^ddy (vertical axis) as a function 
of the ratio E osc /E tot& \ (horizontal axis) for three sets (represented by three 
curves) of sample flows (1) for oj — 1. Dots show computed values of magnetic 
eddy diffusivity. Two dots marked by arrows have been obtained for sample 
flows, for which results of computations are shown on Fig. 3. 

Though time periodicity may be expected to enhance chaotic properties of 
flows (which are necessary for fast dynamo action), our results suggest that 
time dependence such as (1) is not favourable for generation of magnetic 
field. However, dynamos considered here are slow, and therefore this does 
not represent a formal contradiction. 

c) How does magnetic eddy diffusivity depend on temporal frequency of the 
flow? Computations have been carried out for two flows (Fig. 3), which have 
the same profiles of the constituent fields U(x), V c (x) and V s (x). Two curves 
on Fig. 3, computed for the two flows, are labelled (a) and (b) in agreement 
with labelling of the two respective points on Fig. 2 (indicated on Fig. 2 
by arrows). Different patterns of behaviour of magnetic eddy diffusivity are 
observed, though the two flows differ only in amplitudes of their steady and 
time-dependent parts. When uj — > oo, magnetic eddy diffusivity approaches 
finite limit values which can be significantly higher than min w ?7 eddy (u;), as it 
happens in the case (b) (in this case the sample flow has a relatively larger 
time-dependent part of the flow). The plots show that the minimal value of 
7/eddy is achieved at moderate values of temporal frequency, uj = 0(1). 
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Figure 3. Minimal magnetic eddy diffusivity, 7/eddy, (vertical axis) as a func- 
tion of temporal frequency of the flow, u, (horizontal axis). Dots show com- 
puted values of magnetic eddy diffusivity. Curves are labelled (a) and (b), 
referring to the labels of the respective points for uo = 1 on Fig. 2. 



An example of a magnetic mode structure is shown on Fig. 4 and 5. The 
leading term in the expansion of a magnetic mode in the ratio of spatial 
scales, e, is 

Ho = e^ x ^h + ^h fc S fe (x,t)^ 

(see (A6), (A12) and (A22) in the Appendix), where Sfc(x, t) is the solution 
to the first auxiliary problem (11), and h is an eigenvector of the problem 
(13)- (14). The fluctuating part of the magnetic mode, 

G = ]Tn fc S fe (x,t), 
k=i 

is shown on Figs. 4 and 5 as surfaces of constant magnetic energy |G| 2 . 
The Figures show the magnetic mode generated by that flow, for which the 
curve (b) is shown on Fig. 3, for u = 1.5, where the minimum of magnetic 
eddy diffusivity is located on this curve. Cigar-like magnetic structures, 
seen on the Figures, can be associated with stagnation points of the flow 
(1), though no flux-rope solutions similar to those of Zheligovsky (1993b) 
and Galloway and Zheligovsky (1994) are available for time-dependent flows. 
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Figure 4. Isosurfaces of magnetic energy |G| 2 of the fluctuating part of a 
magnetic mode at the level of 40% of the maximal energy. Snapshots of a 
flow periodicity cube (drawn in solid lines) step T/6 are presented. The lower 
vertex is at the point (— ir/2, — tt/2, —tt/2). Dashed lines show the elementary 
cube of the flow stagnation points mesh (toi7T, m 2 7r, m 3 n). 
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Figure 5. Same as Fig. 4, but isosurfaces of magnetic energy of the fluctuating 
part of the magnetic mode at the level of 75% of the maximal energy are 
shown. 
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(Note, that due to parity invariance (2) and space periodicity of the flow 
with the period 2n along each Cartesian axis, the flow velocity vanishes at 
the points (toi7t, 171211, m^Ti), where mi,m2,m3 are arbitrary integers.) The 
sharpness and the shape of the cigar structures vary considerably within the 
time period. 

3. Magnetic eddy diffusivity of flows (1) in the high frequency limit 

Figure 3 suggests, that in the limit uj — > 00 a contribution to magnetic 
eddy diffusivity from the small-scale dynamics of a flow (1) is finite. The 
limit of high frequencies is studied in this Section. 

It is convenient to express the flow (1) in the form 

v(x, t) = U(x) + yfr (V(x)e^ + V(x)e-^*) , (17) 

where V = (V c + iV s )/2. Denote by £ the parabolic magnetic induction 
operator, acting in the space of 27r/u;-periodic small-scale fields F(x, t) with 
a zero mean ((F) = 0): 

dF 

CF = -— + r?V 2 F + V x (v x F). 

ot 

Consider Fourier series of CF: 

00 

CF = CjFe iiut . 

j=-oo 

For the flow (17) 

CjF = -ijuFj+rj^Fj+VxCUxF^+y/u (V x F H + Vx F J+1 ) , (18) 

where Fj are temporal Fourier components of F: 

00 

F(x,f)= £ F,(x)e^. 

j=-oo 

In terms of the Fourier components, the first auxiliary problem takes the 
form 

£jSk = (^U + yfc (s{V + 6Liv)) , (19) 

where S is the Kronecker symbol, and the second one - 

CjT^ = -2r]-^ + U m (S jjk + fe) - US™, (20) 

(V m (S,--!,* + 5{e k ) + r (S j+1> , + - VS£_ 1>fc - VS™ 1>fc ) , 

where Sj^ are Fourier components of S^: 

00 
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Assuming a power-law asymptotic behaviour of the Fourier components 
of Sk and T mjk , one finds 

to be consistent with (18) and (19), and 

r^^^^i/^ + o^-d^ 1 )/ 2 ) 

is consistent with (18) and (20). Here Tj, m ,k are Fourier components of r m fc : 

oo 

r m , fc (x,*)= E r,,m,k(x)J iut - 

j=-oo 

Leading terms of the Fourier series of S k and T myk satisfy 
r]V 2 s 0)k + V x x (2Re(V x si >fc ) + U x s , fe ) = 



si,* = -i + V x x (V x so, fc ) j ; 



(21) 



^V 2 7 , m , fc + V x x (2Re(V x 7l>ro>fc ) + U x 70iro>fc ) 

<9s c 

dx 



-2^ + 2Re(v m s lifc -Vs-), (22) 



li,m,k = -i (V m (s , fc + e fc ) - Vs™ fc + V x x (V x 7 o, ro , fc )) . 
Hence, for uj — > oo coefficients of the eddy diffusivity tensor are 

D m>fe = 2Re (V x 7ljm>fc ) + (U x 70jm>fc ) + 0{uj- 1/2 ). 



The limit values of ~D m ^ are not affected by phase shifts: (21) and (22) 
imply that if V is modified to become e iCt V, where a is a real constant, then 
s ±1)fc changes to e ±ia s±i )fe , 7±i, m , fe - to e ±ia ^ ± ^ m ^ and s 0)fc , 7 0jm)fc and & mik 
remain unaltered. s ^ k and 7o, m ,fc are real; if also V is real (i.e. if V s = 0), 
then s ±l fc and 7±i, m ,jt are imaginary, and hence 

Re(V x si, fe ) = Re(V x 7l)TO)fc ) = 0. 

Therefore, if ImV = 0, then in the limit uj — > oo a contribution from the 
time-periodic part of the flow (17) vanishes. Combined together, these two 
observations imply that a non-zero contribution from the time-periodic part 
of (17) requires linear independence of vector fields V c and V s . 

The same value of molecular viscosity, r\ = 0.1, has been employed in sim- 
ulations. Short-scale magnetic modes can be analyzed in the limit uj — > oo 
along the same lines, as above: their temporal Fourier components exhibit 
the same power-law asymptotical behaviour, as solutions to the auxiliary 
problems. The limit small-scale magnetic modes turn out to be eigenfunc- 
tions of the linear operator defined by the left-hand sides of (21) and (22). 
It has been verified that for each sample flow (1), for which a limit magnetic 
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eddy diffusivity value is reported here, the limit magnetic induction operator 
acting in the space of 27r-periodic small-scale zero-mean magnetic fields has 
no eigenvalues with a positive real part, i.e. for the employed flows rj = 0.1 is 
above the threshold for the onset of generation of small-scale magnetic fields 
in the limit of high temporal frequencies. 

We carried out the following numerical experiments: 

a) Limit values of magnetic eddy diffusivities have been computed using 
(21) and (22) for those three sets of sample flows (1) (see Fig. 6), for which 
plots of Fig. 2 have been constructed for to — 1. The left-most points (for 
-E'osc — 0) of two respective curves on Figs. 2 and 6, obtained for the same set 
of sample flows, show the same value; this can be used to identify pairs of the 
respective curves. Like for finite frequency (see Section 2), a relative increase 
of energy contained in the time-dependent part of the flow is accompanied 
by an overall increase of magnetic eddy diffusivity, though the dependence is 
not necessarily monotonous. However, the influence of time-dependent parts 
of flows has now weakened. 




0.0 0.2 0.4 0.6 0.8 1.0 

Figure 6. Minimal magnetic eddy diffusivity ^ c dd y (vertical axis) in the limit 
uj — > oo as a function of the ratio E osc /E tota i (horizontal axis) for three sets 
(represented by three curves) of sample flows (1). Dots show computed values 
of magnetic eddy diffusivity. 
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b) Distribution of minimal magnetic eddy diffusivity for (1) has been exam- 
ined under the conditions, apparently least favourable for generation - in the 
absence of the steady component (i.e. for U = 0) and in the limit to — > oo. 
A histogram of limit values of ?7eddy computed for 45 sample flows, satisfying 
(15), (16) and U = 0, is shown on Fig. 7. Only in 2 cases out of 45 the limit 
magnetic eddy diffusivity is negative. 
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Figure 7. Histogram of minimal magnetic eddy diffusivity values in the limit 
uj — > oo for 45 sample flows satisfying (15), (16) and U = 0. 

Conclusion 



Generation of magnetic field involving large spatial scales by time- and 
space-periodic small-scale parity-invariant flows has been studied. The aniso- 
tropic magnetic eddy diffusivity tensor has been calculated. A complete ex- 
pansion of magnetic modes and their growth rates in power series in the scale 
ratio has been constructed for flows (1). Simulations have been carried out 
for flows (1) with random harmonic composition and an exponentially decay- 
ing energy spectrum. Flows giving rise to negative magnetic eddy diffusivity 
have been found, for molecular diffusivity above the instability threshold for 
small-scale magnetic field generation. Numerical results suggest that gener- 
ation of large-scale magnetic field by time-periodic flows is less efficient than 
generation by steady flows of the same energy. High temporal frequencies 
may also be unfavourable for generation. Nevertheless, we have found nu- 
merically instances of flows (1) with a vanishing steady component, for which 
magnetic eddy diffusivity is negative for u — > oo. 
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Appendix. Formal asymptotic decomposition of magnetic modes 
and of their growth rates for time-periodic flows (1) 

We derive complete asymptotic expansions of a magnetic mode and its 
growth rate for the flow (1), which it is convenient to express here as 

v(x, t) = U(x) + W(x)e^* + W(x)e _iwt , (Al) 

where the notation 

W = v^JV (A2) 

is used. Our derivation follows that of Zheligovsky et al. (2001) for steady 
flows. 

As in the main body of the paper, x G R 3 denotes the fast, and y = ex 
- the slow spatial variable, (■) and {■} denote the mean and the fluctuating 
part of a vector field, respectively: 

(f(x,y,t)) = (27r)- 3 / f(x,y,t)rfx, {f (x, y, t)} = f - (f) . 

The following assumptions are made concerning the flow v: U and W 
are 27r-periodic in fast Cartesian variables, independent of time and of the 
slow variables, solenoidal: 

V-U = V-W = 0, (A3) 

and parity-invariant: 

U(x) = -U(-x), W(x) = -W(-x). (A4) 

A magnetic mode H(x, y,t) is a solenoidal (5) solution to the Floquet 
problem (4), which is 27r-periodic in each spatial variable and has the same 
temporal period T = 2ir/uj, as the flow. For the flow (Al), Fourier compo- 
nents of a magnetic mode 

oo 

H(x,y,t)= ]T H,(x,y)e^ 

j=-oo 

satisfy 

XHj = -ijullj + r^V 2 ^- + V x (W x lUxH^Wx H J+1 ) (A5) 
for all j. 

A solution to this system of equations is sought in the form of power 
series 

oo 

H,(x,y) = E(%n(y) + G,>(x,y))e n (A6) 

n=0 

and (7). In (A6) Hj >n an d Gj, n are the mean and the fluctuating part of the 
respective term of the series: (Gj, n ) — 0. Obviously, any term H n of the 
expansion (8) can be recovered from (A6): 



H n (x,y,t)= J2 (H j , n (y) + G J> (x,y))e^ 

j=-oo 
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After modification of the spatial gradient (6), expansion and separation 
of mean and fluctuating parts of each Fourier harmonics, the solenoidality 
condition (5) reduces to 

V y • H jn = 0, (A7) 
V x • Gj, n + V y • Gj, n -1 — 

for all j and n > 0, where it is understood Gj, n = for n < 0. Here and in 
what follows the subscripts x and y refer to differential operators in fast and 
slow variables, respectively. 

Substitution of (A6), (7) and (6) into (A5) yields 



E 

n=0L 



ijuH jin + CjG. tn + rj (2(V X • V y )G J - n _i + V y (H jjn _ 2 + G,- n _ 2 )) 



+V X x (w x Hj_i >n + U x H j>n + W x H j+1 
+V y x (W x (Hj_i, n _i + Gj_i, n _i) + U x (H j>n _i + Gj-n-i) 

n 

+W x (Hj + i jn _i + Gj + i )Tl _i)) — A n _ m (H_ 7 - )jn + Gj jTn ) 



m=0 

Here it is denoted 



e n = 0. (A8) 



£jF = 77V x Fj + V x x (W x Fj_i + UxF 3 +Wx F i+1 ) - ijuFj, 
where Fj are Fourier components of 

oo 

F(x,y,t)= £ F^x.y)^. 

j=-oo 

We make a final technical assumption that for any function f(x, t), 
T-periodic in time and 27r-periodic in space, such that (f) = for all t, 
the problem CF = f has a unique small-scale solenoidal solution with a van- 
ishing spatial mean, which has the same time and space periodicities, as the 
flow. ( Equivalent ly, the magnetic induction operator £ is assumed to have a 
trivial kernel.) Generically this condition holds. 

We proceed by successively equating the mean and the fluctuating part 
of each term of the series (A8) to zero. 

%. The leading (n = 0) term of (A8) takes the form 

CjG.,0 + (H,-_ li0 • V X )W + (H,- o • V X )U + (H j+1>0 • V X )W 

= A (H ii0 + Gj-o) + ijuHjfl. (A9) 

The mean of (A9) is 

= (iju + Ao)H„ . (A10) 

Thus it can be assumed 

H il0 = VjVO; A = (All) 
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(this is a normalisation condition: any other formal solution to (A10): 
Hj t o = Vj 7^ J, A = — iJoj for J ^ 0, does not represent any new eigen- 
solution to the original Floquet problem, since solutions H to the Floquet 
problem (4) are defined up to a factor Ce lJu)t ). 

Consequently, the fluctuating part of (A9) yields by linearity 

G,,o = ES^(x)H fc (y), (A12) 
k=i 

where vector fields S^* satisfy 

CjS., k = -JL (5{W + 5iU + Si.W) , (A13) 

which is a representation of the first auxiliary problem (19) involving notation 
(A2). The problem (A13) has a unique solution by the assumption that the 
kernel of C is empty. It is evident from (A13) and from the definition of the 
operators Cj that 

S -j,fc = S i,fc; V x • S jtk = Vj, k. 

Parity invariance of the flow (2) implies that parity anti-invariant vector 
fields are an invariant subspace of C, and since the right-hand side of (A13) 
is parity anti-invariant, so are S^: 

S^(x) = S i>Jfc (-x). (A14) 

Divergence (in fast variables) of (A13) implies that are solenoidal. 

%%. The second (n — 1) term of (A8) reduces with the use of (A3), (A7) and 
(All) to 

CjG. tl + 2^(V X • V y )G ji o + (H.j-1,1 ■ V X )W + (H jA ■ V X )U + (H j+M • V X )W 
+V y x (W x Gj-1,0 + U x Gj, o + W x G J+1)0 ) 

- ((^W + 5iU + (JiiW) • V y ) H ,o = A^^Hco + G jt0 ) + iju>H jA . (A15) 
Upon substitution of (A12) the mean of (A15) becomes 

3 

V y x^(Wx Sj_i )fc + U x Sj, fc + W x S j+ i ifc ) Hg j0 = Ai^H ,o + ijwH^. 
fc=i 

In view of (A4) and (A14), the averaged cross products at the left-hand side 
of this equation vanish and thus 

H jil = VjVO; A 1= 0. (A16) 

After (A16) and the representations (A12) are plugged in, the fluctuating 
part of (A15) becomes 
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3 3 

+ EE 

k=l m=l 



£,0,! = -(H 0) i • V x ) (<S'W + $U + SL^W) 
- 2 V -^ + W m (s,-_ 1>fc + fe) - WS£ 1>fc 



(A17) 



+U m (S jife + <Sfc fc ) - US™ + F (S j+1>fc + - WS" 

Hence by linearity 



m 
+l,fc 



dy r , 



3 

E 

fc=l 



3 3 

EE 

fc=l m=l 



G jtl = £ S ilfc (x)H* x (y) + £ E IW(x) — M( y)> (A18) 



where vector fields I\,- jmjfc (x) satisfy 



dS 



CjT.^ = -Zv-gf 1 + W m (S,_ 1)fe + 5{e k ) - WSf_ ljfc 



+U m (Sj, fc + fe) - US™ fc + W™ (S j+1)fe + - WS™ ljfc , (Af 9) 

which is a representation of the second auxiliary problem (20) involving no- 
tation (A2). By standard arguments it is verified that 



k ~Fj,m,ki ^j,m,k{p^) ^j,m,k{, x). 



Divergence (in fast variables) of (A19) implies V x • ^j,m,k + S™ fc = 0. 

Hi. The third (n = 2) term of (A8) by virtue of (All) and (A16) reduces to 

£,■0,2 + v (2(V X • V y )G ifl + VjOSgHo.0 + G,- )) 

+(H j _ li2 • V X )W + (H ji2 • V X )U + (H j+li2 • V X )W 
+V y x (W x G,_ M + U x G jt i + W x G j+ i,i) 

- ((5{W + 5iU + 6LiW) • V y ) H 01 = A 2 ($H 0)0 + G jt0 ) + iju>H jt2 . (A20) 
In view of (A18), (A4) and (A14) the mean of this equations is 



3 3 



77^v y H , +v y x E E ( w x r i-i,m,fc + u x r jim>fc + w x r j+ i, m)fc 



m=l fe=l 



<9H 



0,0 



A 2 <^H ,o + ijuH jt 2- 



dy m 
(A21) 



Thus the leading terms of the expansions of the mean magnetic field, H , 
and of the growth rate, A 2 , are a solution to the eigenvalue problem 



3 3 

A4H 0i0 = ^V 2 H ,o + V y x ]T £ (2Re ( W x r 1)m , fc ) + (U x r , m , fc ) 

m=l k=l 



— A 2 Ho,o, 



dy m 
(A22) 
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arising from (A21) for j = 0. (A22) is equivalent to (9) in view of (A2). 
Subsequently one obtains from (A21) 

1 3 3 dH k 

Hj,2 = — V y x^^(Wx Tj-i,™,k + U x Tj^ k + W x T j+ljmi k 



0,0 

for any j 7^ 0. Hence now the quantities Gj are also entirely determined by 
(A12). " ' 

Bounded solutions to the eigenvalue problem (A22) are Fourier harmonics 
H 0j0 = he 4qy , satisfying the orthogonality condition (14) and 

3 

??|q| 2 h + qx^ (2Re (W x r_ 1)TOifc ) + (U x r 0)mjfc ))q m h fc = -A 2 h, (A23) 

m=l 

which is an equivalent of (13) in new notation (A2). 

The x-dependent prefactors in front of unknown vector fields Hq 2 and 
dHQ^/dy m in the fluctuating part of (A20), 

3 ' 2 " ~ v \ hh U^^T + h dx> dy mdyi ) + VyG -°J 

— (H 3 -_i j2 • V X )W - (H j)2 • V X )U - (H j+1>2 • V X )W 

+ ({5{W + 5iU + SUW) • V y ) H ,i (A24) 

" V y x ( £ (W x S,--!,* + U x S jife + W x S j+1 , k ) H fc , 
\fe=i 

+ EE{ Wx r i-i,m,fc + U x T hmtk + W x r i+1 , mife } — -M + A 2 G i;0 , 

are the same as those in front of Hq 1 and dH.Q Q /dy m , respectively, in (A17), 
which by linearity implies a representation 

3 33 <9H fc 

Gi,2 = E S^(x)H fc 2 (y) + EE r,, m , fe (x)— M (y) + Q, i2 (x,y). 
fc=i fc=i m=i a y m 

Vector fields Q J)2 can be found from a system equations, obtained from (A24) 
by changing G. i2 — > Q. 2 and dropping all terms involving 2 or derivatives 
of Hq x (the right-hand sides of the resultant equations are at this stage 
known) . 

iv. Let A' 2 and h' be the second eigenvalue and the associated eigenvector, 
satisfying (A23) and (14). Subsequent (n > 2) terms of (A8) provide a 
hierarchy of equations, which can be solved under the condition A 2 7^ A' 2 . 
Equations for n < N yield: 

• vector fields Hj, n f° r ah j 7^ and n < N; 

• vector fields H ,„ for all n < N — 2; 

• vector fields Gj :Tl for all j and n < N — 2; 
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representations of Gj, n of the form 



3 s 3 dH k 



G jtn = E^W H Uy) + E E r,, m , fe (x)— ^(y) + Q J> (x,y) (A25) 

fc=1 fc=1 m=1 y i/m 



for n = N — 1 and n = N — 2 with known vector fields Qj, n , (Qj>) = 0; 
• quantities X n for all n < N. 

Upon substitution of (A25) for n = N — 1 the mean of the equation 
corresponding to n = N in (A8) becomes 

3 3 ^ QYl k 

^v 2 h,, jv-2+VyX E E ( w x r i-i,m,* + u x r i>mifc + w x r j+hmik ) °> N - 2 
m =i k=i °y m 

+V y x x Qj_i,jv-i + U x Qj-at-1 + W x Q j+1)Ar -i^ 

JV-2 

= ^] AAr_ m H jim + ijuH jjN . (A26) 

m=0 

Consider (A26) for j = 0: 

(M. — A 2 )H 0i Ar-2 — AatH 0i o 

JV-3 

= Ajv-mH 0im - V y x (2Re (W x Q 1)iV -i) + (U x Q ,jv-i>) , (A27) 

m=l 

where the right-hand side is a known vector field. Projecting this equation out 
in the direction of Ho.o one can uniquely determine Xn- In the complementary 
invariant subspace of M. the operator M. — A 2 is invertible, and thus H ,Ar-2 
can be determined from (A27) up to an arbitrary multiple of H 0i o, which we 
can demand to vanish. Now Hj,N for j ^ can be found from (A26), and 
Gj,N-2 are determined by (A25) for n = N — 2. 

The fluctuating part of the equation corresponding to n = N in (A8) 
becomes after the substitution of (A25) for n = N — 1 

3 ' ,N V feTi m^i \ dx m dy m fr{ dxi dyidy m ) 

+2(V X • V y )Qj i AT_i + V*G jiJV _ 2 ) - (H.j-i,N ■ V X )W + (W • Vy)!!^^ 
-(H i>JV • V X )U + (U • V y )H,. N _, - (H i+liJV ■ V X )W + (W ■ V^H,-^-! 

-Vy X 



OJV-2 



E( Wx S,-i, fe + U x S hk + W x S i+1>fc ) H fe N _, (A28) 
_fc=i 

+ E E { w x ri-i.rn.fc + u x r,- m , k + w x r i+1>m>fc } 

fc=i m=i 

1 iV-2 

+ {W x Qj-i.AT-i + U x Q jiJV -i + W x Q j+ i,jv_i)} 



m=0 
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Like in the case of (A24), the structure of this equation implies by linearity 
the representation (A25) of G^n, where the quantities Qj,jv are uniquely 
defined by 

/ 3 3 3 o-p /)2TTfc N 

c jQ ., N = -J 2 e E E £*\,°l- 2 + 2 ^ ■ V -)Q— + v y G — 



k=l m 



-(1 - (Hj_ 1)JV • V X )W - (1 - <Jg)(H^ • V X )U - (1 - 5i 1 )(H 3 - +lj jv • V X )W 
+(l-5D(W-V y )H J _ 1 ^_ 1 +(l-^XU-V y )H^_ 1 +(l-5iO(W-V y )H J+liAr _ 1 



-V y x 



3 3 Qftk 

E E { w >< r J-i,m,k + u x r i>m>fc + w x r i+1>m>fc } 

.k=i m=i 



+ {W x Qj_i >JV _i + U x Qj- jv-i + W x Q j+1;Ar _i)} 



JV-2 
m=0 



(the right-hand side of this equation is known). This equation was obtained 
by omitting in (A28) all terms, involving H 0i at or derivatives of H 0; at-i, and 
changing G. jN to Q.at. 

Thus a complete asymptotic expansion of magnetic modes and their 
growth rates is constructed. Like in the case of a stationary velocity, it 
can be easily verified that 

H 0i0 = he** H ,„ = Xnh'e 1 ^ Wn > 0, H jn = h^e™ Vn, j + 

(where h Jjra =const, hj j7l • q = 0), 

G Jin = g„(x)e^ Vn>0 

and thus the eigenmode admits a representation 

H = e ieqx h(x, t), q = const. (A29) 

(For this reason for constructions of section iv it was sufficient to demand 
that A' 7^ A2, and not that M. — A2 is invertible in the whole domain.) This 
stems from the fact that for the velocity (Al) the domain of the magnetic 
induction operator splits into invariant subspaces, each comprised of vector 
fields (A29) and categorised by wavevectors q. 

Analyzing parity of solutions of the hierarchy of equations (A8) con- 
structed in this Appendix, one finds that all Gj, n with even indices n are par- 
ity anti-invariant; all Gj jTl with odd indices n are parity- invariant; (Hj j7l ) = 
for any odd n; and A n = for any odd n. 
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